Pulsar kicks from a dark-matter sterile neutrino 
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We show that a sterile neutrino with mass in the 1-20 keV range and a small mixing with the 
electron neutrino can simultaneously explain the origin of the pulsar motions and the dark matter 
in the universe. An asymmetric neutrino emission from a hot nascent neutron star can be the 
explanation of the observed pulsar velocities. In addition to the pulsar kick mechanism based on 
resonant neutrino transitions, we point out a new possibility: an asymmetric off-resonant emission 
of sterile neutrinos. The two cases correspond to different values of the masses and mixing angles. 
In both cases we identify the ranges of parameters consistent with the pulsar kick, as well as 
cosmological constraints. 



Pulsars are known to have large velocities ranging from 
100 to 1600 km s _1 0. In contrast, the average velocity 
of an ordinary star in our galaxy is 30 km s . Born in a 
supernova explosion of an ordinary star, the pulsar must, 
therefore, receive a substantial "kick" at birth. The high 
angular velocities of pulsars may also be explained by the 
same kick 2] ■ The origin of the pulsars kicks remains an 
intriguing outstanding puzzle. 

According to numerical simulations of supernova ex- 
plosions |3j, the asymmetries in the core collapse could 
not account for a kick velocity of more than 300- 
600 km s" 1 . Although the average pulsar velocity is in 
this range |l| , there is a substantial population of pulsars 
with velocities in excess of 700 km/s, including some with 
speeds over 1000 km s . These velocities appear to be 
too large to arise from an asymmetry in convection. Evo- 
lution of close binary systems has also been considered 
as a source of a kick velocity Q. Alternatively, it was 
argued [j| that the pulsar may be accelerated during the 
first few months after the supernova explosion by its elec- 
tromagnetic radiation, the asymmetry resulting from the 
magnetic dipole moment being inclined to the rotation 
axis and offset from the center of the star. Both of these 
mechanisms, however, have difficulties explaining the ob- 
served magnitudes of the pulsar kick velocities. 

In this paper we point out a new way in which the 
pulsar kicks could be generated from off-resonant active- 
sterile neutrino conversion in the neutron star core. In 
addition, we identify the parameters for which the reso- 
nant conversions can take place. Both cases are consis- 
tent with the sterile neutrino being the dark matter. 

The neutrinos emitted from a nascent hot neutron star 
carry a total energy of 10 53 erg. Only a 1% asymmetry in 
the distribution of these neutrinos could produce a kick 
consistent with the observed magnitudes of pulsar ve- 
locities. Despite negligible magnetic moments, neutrinos 
are affected by the strong magnetic field of the neutron 
star through their interactions with the polarized elec- 
trons in the medium. In particular, resonant neutrino 
oscillations occur at different densities depending on the 



direction of the magnetic field relative to the neutrino 
momentum 0, Q ■ Oscillations of active neutrinos could 
explain the observed velocities of pulsars if the resonant 
conversion f MjT <-> v e to ok p lace between two different 
neutrinospheres [^.l9lllollllLlT^|. However, the neutrino 
masses required for the resonant transition between the 
two neutrinospheres, at density 10 11 — 10 12 gcm -3 , are 
inconsistent with the present limits on the masses of stan- 
dard electroweak neutrinos. 

These limits do not apply, however, to sterile neutri- 
nos that may have only a small mixing with the ordinary 
neutrinos. Resonant oscillations of active neutrinos into 
sterile neutrinos have been pro posed as a possible expla- 
nation of the pulsar kicks [lj| ■ It was recently pointed 
out that conversions of trapped neutrinos in the core of a 
hot neutron star can make the effective mixing angle close 
to that in vacuum pdj ]. This has important implications 
for the pulsar kick mechanism because, in the absence 
of matter suppression of the mixing, the non-resonant 
production of sterile neutrinos can dominate. This is 
asymmetric, and it can explain the pulsar velocities as 
well. The same sterile neutrino in the keV mass range 
can be the dark matter in the universe [IdL llj| . Ster- 
ile neutrino with the requisite masses and mixing appear 
in a number of theoretical models |16l. for example, in 
models with broken mirror parity [171 ]. 

We will assume that neutrinos have negligible magnetic 
moments and that they have only the standard interac- 
tions with matter. Unusually large magnetic moments of 
neutrinos 0], as well as other kinds of new physics 
have also been considered as the origin of the pulsar kicks. 
We will further assume that only v e has a significant mix- 
ing with u s , characterized by sin 2 9 ~ 10~ n — 10~ 7 , while 
the oscillations between v^ T and v s are suppressed by 
small mixing angle and/or matter effects. 

For a sufficiently small mixing angle between v e and v s , 
only one of the two mass eigenstates, v\, is trapped. The 
orthogonal state, \v-i) = cos0 m \v s ) +sin0 m |i/ e ), escapes 
from the star freely. This state is produced in the same 
basic urea reactions (i> e + n ^ p+e~ and i>e+P ^ n+e + ) 
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with the effective Lagrangian coupling equal the weak 
coupling times sin 9 m . 

Urea processes are affected by the magnetic field, and 
so the active neutrinos are produced asymmetrically de- 
pending on the direction of their momenta relative to the 
magnetic field. Chugai [2(j and Dorofeev et al. [2l| have 
proposed that this asymmetry might explain the pulsar 
kick velocities. However, the asymmetry in the produc- 
tion amplitudes does not lead to an appreciable asymme- 
try in the emission of neutrinos because this asymmetry 
is washed out by scattering 123 ■ Therefore, the mecha- 
nism considered in Refs. [20L l2l| does not work. 

However, since 1 1^2} is not trapped anywhere in the 
star, the asymmetry in its emission is exactly equal to 
the production asymmetry. Depending on the parame- 
ters, this asymmetry can be as large as 25%, just as in 
the case of the active neutrinos |2j . If the sterile neu- 
trinos carry away about 10% of the thermal energy, the 
resulting overall asymmetry can reach the required few 
percent, which is what one needs to explain the pulsar 
kick velocities. 

The production cross section of v% depends on the ef- 
fective mixing angle in matter m , which, in general, is 
different from the vacuum mixing angle 0: 



sin 2 20™ = 



(Am 2 /2p) 2 sin 2 20 



(Am 2 /2p) 2 sin 2 20 + (Am 2 /2p cos 20 - V m ) 2 ' 

where the matter potential V m is positive (negative) for 
v(y)i respectively; p is the momentum. For the case of 
v e , which is relevant for our discussion, V m reads: 
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In a core collapse supernova, the initial value of this 
matter potential is V m ~ (—0.2... + 0.5)Vb, where Vb = 
G F p/V2m n ~ 3.8 eV(p/10 14 gcm- 3 ). The average prob- 
ability of v e — ► v s conversion in presence of matter is 



(Pm) = \ 
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where we have included a possible suppression due to 
quantum damping 23] , which depends on the oscillation 
length A osc and scattering length A s . 

It was pointed out in Ref. [14| that, in the presence 
of sterile neutrinos, rapid conversions can take place be- 
tween different neutrino flavors. This, in turn, can drive 
the effective potential to its stable equilibrium fixed point 
V m — > 0. This equilibration takes place on a very short 
time scale and results in the destruction of the initial 
asymmetry of v e over v e . (The initial electron lepton 
number asymmetry in the supernova core depends on 
the details of the prior stellar collapse history and, es- 
pecially, the electron capture on heavy nuclei and neu- 
trino transport in the in-fall epoch.) This, in turn, leaves 
the system out of /3-equilibrium. On a time scale ~ms, 



electron capture reactions in the core will return the sys- 
tem to beta equilibrium and the v e so produced will be 
converted. These processes will continue until both beta 
equilibrium and a steady state equilibration arc achieved. 
Once the equilibrium is achieved, the effective mixing an- 
gle in matter is close to that in vacuum. In Ref. 0] 1 this 
effect was not taken into account, and only the resonant 
emission from a thin shell was considered. A larger mix- 
ing angle makes a big difference in that the emission off 
resonance becomes possible from the entire volume of the 
neutron star core. 

The equilibration mechanism relies on the fact that, in 
general, the probability of v e — > v s oscillations is different 
for neutrinos and antineutrinos due to the opposite sign 
of V m in eq. QJ [3^J . The resulting change in the amounts 
of v e and v e drives V m to zero. Following Ref. [14J, we 
estimate the time scale for this process: 
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where dli = (2ir 2 )~ 1 e 2 de v , and Vm^ is the initial value 
of the matter potential V m . 

Since the neutrino emission depends on the value of 
(P m ), the time scales for resonant and off-resonant con- 
versions differ. In the vacuum case, Am 2 /(2(E)) > \V m \, 
this mechanism is not relevant, because V m is negligible 
from the start. In the case of resonance, Am 2 /(2(E)) ~ 
V m , the conversions of v^s are strongly enhanced over 
the Ve§- Therefore the equilibration mechanism is very 
efficient. Using eq. |@}, we estimate the time: 
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This holds as long as the resonance condition V r 
Am 2 1 (2(E)) is satisfied. The back-reaction effects on 
V m due to the change in Y Ve may slow down the equili- 
bration. Their detailed analysis is beyond the scope of 
this paper; we do not expect the back-reaction to alter 
the time scales significantly. 

If neutrinos are converted off-resonance, it is necessary 
to evaluate the integral in eq. (@J taking into account 
that both the chemical potential and the reaction rates 
for neutrinos and antineutrinos are different. An approx- 
imate evaluation of the time scale yields: 
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FIG. 1: Neutrino conversions can explain the pulsar kicks for 
values of parameters either in region 1 (Ref. Ill), or in region 
2 [present work]. In a band near the f2„ = 0.3 line, the sterile 
neutrino can be the dark matter. 



Requiring that the off-resonance equilibration takes place 
during the first second, imposes a combined lower bound 
on Am and sin 29. In Fig. □ this bound, for = 
0.1 eV, fi = 50 MeV, T — 20 MeV, limits region 2 from 
below. (T = 20 MeV is plausible for the immediate post 
core bounce conditions, where models suggest a core tem- 
perature between tens of MeV and about 60 MeV |24|-) 

The probability of oscillation depends on the values 
of Am 2 , sin (9, and V m . One can identify three different 
regimes, (i) the "vacuum" case: Am 2 /(2(E)) \V m \; 
(ii) the resonance regime, Am 2 / (2(E)) w V m , in which 
the conversion probability is strongly enhanced; (iii) the 
"matter-suppressed" oscillations, for Am 2 /(2(E)) <C 
V m , which have a vanishing probability. 

Whether or not a resonance occurs anywhere in the 
core depends on the mass of the sterile neutrino and its 
mixing angle with v e . The resonance occurs at some 
point in the core for parameters in the range marked "1" 
in Fig. ^ The lower boundary of this region comes from 
enforcing the adiabaticity and the weak damping condi- 
tions. The vertical span reflects the range of densities in 
the outer core of a neutron star. 

The case of resonant emission. If the resonance 
condition is satisfied (region 1 in Fig. ^) , the emissivity 
is strongly enhanced in the region of resonant density, 
while it is suppressed elsewhere. The resonant region is a 
thin layer, whose deviation from a spherical shapejs_due 
to the electron polarization in the magnetic field 
IT^ |. For small mixing angles (region 1), the time r° ff " 
is longer than a few seconds for some reasonable initial 
conditions. Hence, the matter potential is not affected 
outside the resonant region. Therefore, in region 1 the 
pulsar kick mechanism works as described in Ref. |13| . 
where the equilibration was not considered |33|. 

Off-resonance emission after V m has equili- 
brated to zero. To get the pulsar kick, we require that 



a substantial fraction of energy be emitted in sterile neu- 
trinos. The dominant process of neutrino production in 
the core is urea. One can estimate the emissivity in ster- 
ile neutrinos following Ref. (l4| : 
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The emissivity depends on the value of (P m ). In the 
"vacuum" regime the emissivity is proportional to sin 2 26 
and is independent of Am 2 . The total energy emitted per 
unit mass is simply the emissivity multiplied by the time 
of emission, which is 5-10 s. For given values of \i and 
T , the requirement of emitting a certain fraction of the 
energy in sterile neutrinos translates into a lower bound 
on the allowed values of sin 2 29, as shown in Fig.Q] Such 
a bound is independent of Am 2 , as long as the V m ~ 
condition is satisfied. If the equilibration process is effi- 
cient, V m w and the matter effects are negligible. Then 
the oscillations v e — > v s are controlled by the vacuum 
mixing angle sin 2 29. 

As discussed by Dorofeev et al. [2lJ i the urea processes 
result in an asymmetric neutrino production if the elec- 
trons are polarized. This asymmetry arises from the spin- 
multiplicity factor. Only the electrons in the lowest Lan- 
dau level (n = 0) contribute to this asymmetry. The 
number of neutrinos dN emitted into a solid angle dfl 
can be written as ^ = N (l + ecosO„), where 0„ is 
the angle between the direction of the magnetic field and 
the neutrino momentum, and Nq is some normalization 
factor. The asymmetry parameter e is equal 
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where g v and g A are the axial and vector couplings, £t Q t 
and £ s are the total neutrino luminosity and that in ster- 
ile neutrinos, respectively. Following Ref. |^J, we obtain 
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Here n is the electron Landau level, and we have as- 
sumed that protons and neutrons are non-relativistic, 
non-degenerate, and that their polarization can be ne- 
glected. The effective electron mass at high density, mt , 
remains below 10-11 MeV everywhere in the star |24tl25| . 
Its exact value does not affect the expression for F(n) in 
ea. dlOJl by more than a few per cent. If only the lowest 
Landau level is populated, then ko = 1. The weaker the 
magnetic field, the more levels are populated, the smaller 
is the fraction fc . In Fig. [21 we show ko as a function of 
chemical potential y, for different values of the magnetic 
field B. We took the temperature to be T = 20 MeV. 
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FIG. 2: The fraction of electrons in the lowest Landau level as 
a function of fi for T — 20 MeV. The value of the magnetic 
field in the core of a neutron star is shown next to each curve. 

The ratio (£ s /£tot) = £ s /{£ s + £u) can be estimated 
from the following: 

(%\~stf6 f^-) 6 / M / d .o,., (11) 

J \ ^ v— sphere / 

where f M < 1 is the fraction of enclosed mass of the 
core from which the emission of sterile neutrinos is effi- 
cient, and /d.o.f. > 1 is an enhancement due to a possi- 
ble increase in the effective degrees of freedom at high 
density j24[. We take the core temperature T corc w 
(20 — 50)MeV, while the neutrinosphere temperature is 
H_ B phere « (2 - 10)MeV. Hence, for sin 2 6 ~ 10~ 8 the 
ratio {£ s /£tot) can be in the range 

r £ = ~ 0.05 - 0.7. (12) 

\<-Hot / 

A higher ratio is in conflict with the observation of neu- 
trinos from SN1987A, while a lower value does not give a 
sufficient pulsar kick. Our point here is that, for reason- 
able values of parameters, the range in eq. I|12H is possible, 
even if not necessary. 

Finally, the magnitude of the asymmetry is 




which can be of the order of the requisite few per cent for 
magnetic fields 10 15 — 10 16 G, as can be seen from Fig. El 

The magnetic field at the surface of an average pul- 
sar is estimated to be of the order of 10 12 — 10 13 G [26j. 
However, the magnetic field inside a neutron star may 
be as high as 10 16 G H EHH HJ. The existence of 
such a strong magnetic field is suggested by the dynam- 
ics of formation of the neutron stars, as well as by the 
stability of the poloidal magnetic field outside the pul- 
sar [23. Moreover, the discovery of soft gamma repeaters 
and their identification as magnetars [28j, i.e., neutron 
stars with surface magnetic fields as large as 10 15 G, 
gives one a strong reason to believe that the interiors of 
many neutron stars may have magnetic fields as large as 
10 15 — 10 16 G and that only in some cases this large mag- 
netic field breaks out to the surface. There is also a plau- 
sible physical mechanism by means of which such a large 
magnetic field can form inside a neutron star through the 
so called a — fl dynamo effect 29, 30j. Nuclear matter 
inside a cooling neutron star is expected to have a suf- 
ficient amount of convection to generate magnetic fields 
as large as 10 16 G 0. 

The ratio of luminosities rg translates into an upper 
bound on the mixing angle. The allowed region for the 
masses and mixing angles is marked "2" in Fig. ^ 

Gravity waves from a pulsar kick due to neutrino con- 
versions might be observable }31j and could provide a way 
to test this mechanism. 

To summarize, we have analyzed the pulsar kick re- 
sulting from active to sterile neutrino conversions. We 
have shown that, in addition to resonant transitions |13| . 
the off-resonant neutrino oscillations can also result in a 
pulsar kick consistent with observations. The two mech- 
anisms require different masses and mixing angles. We 
have identified the range of parameters for which each of 
the two mechanisms can explain the observed velocities 
of pulsars. Part of this range is consistent with the sterile 

WDM 

neutrino being the dark matter with f2„ sa 0.27. 
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